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1. Introduction 

In two Euclidean dimensions the eigenvalue distribution of the SU(Af) Wilson matrix associ- 
ated with a non-selfmtersecting loop undergoes a phase transition in the infinite -N limit as the loop 
is dilated [1]. This phase transition has universal properties shared across dimensions and across 
analog two-dimensional models [2, 3]. Thus, a detailed understanding of the transition region in 
2D is of relevance to crossovers from weakly to strongly interacting regimes in a wide class of 
models based on doubly indexed dynamical variables with symmetry SU(^V). 

We are focusing on the eigenvalues of the Wilson loop. The associated observables are two 
different functions p^' ue (0), p^ m (8) of an angular variable 6. At infinite N the two functions have 
identical limits: p£ ue (0) = p~ ym (0) = 0.(0). 

For a specific critical scale, the nonnegative function p=o(0) exhibits a transition at which a 
gap centered at 6 = ±n, present for small loops, just closes. This transition was discovered by 
Durhuus and Olesen in 1981 [1]. 



2. Eigenvalue densities 

The probability density (with respect to the Haar measure) for the Wilson loop matrix W is 
given by 

&N(W,t) = Y,drXr(W)e-^ c ^ (2.1) 

r 

with t = Xsrf, where A = g 2 N is the standard 't Hooft coupling and srf denotes the area enclosed by 
the loop. The sum over r is over all distinct irreducible representations of SU(A r ) (with dimension 
d r and quadratic Casimir C2(r)). Xr(W) is the character of W in the representation r. 
Density functions are obtained from 

GTiz) = I (Tr^) = i^Oogdctfr- W)>, (2-2) 

GTiz) = -^log(det(z- wy l ) (2.3) 



through (£ = true, sym) 



p t N {B)=2 lim Re e w+e G f N (e w+£ " 
e^0+ 



1 . (2.4) 



Only p l ™ e has a natural interpretation at finite N, it literally is the eigenvalue density. p^ ym is 
determined by the averages of the characters of W in all totally symmetric representations. The 
function is monotonic on each of the segments (—71,0) and (0, n). It turns out that the appropriate 
area variable for p^ ym is not t but T = t{\ - 1 /N) [4]. When p^J m (d,T) is compared to p^ e {d,t), 
the 1 /N correction in t relative to T has to be taken into account. 

The infinite-A^ critical point is at T = 4. For T > 4, p^ m (d,T) approaches poo(d,T) by power 
corrections in l/N [4]. For T < 4, p„(d,T) is zero for |0| > 6 C (T), where < C (T) < % and 
C (4) = K (cf. Fig. 1). In this interval p^ m (d,T) approaches zero by corrections that are exponen- 
tially suppressed in N. 

The true eigenvalue density p l ™ e (d,t) has ,/V peaks (in the interval [— %,%]) and oscillates 
around p^ ym (0,r) (cf. [5]). 
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Figure 1: Plots of p^ m (9,T) for T = 2 (left), T = 5 (right), and N = 3,5, 10,25,50, 100,250 together with 
P-(0,T). 

3. Integral representations 

The density p^" e can be obtained from the expectation value of [5] 

R 



(3.1) 

p=0q=0 

with |v| < 1. Xp (W) (resp. %q(W)) denotes the character of W in a totally antisymmetric (resp. 
symmetric) representation. When we set u = — v + £ and expand to linear order in e, the LHS reads 

1 



R{-v + e,v,W) = 1 -eTr 



(3-2) 



v-W*' 

After decomposing the tensor product p A ® q s into irreducible representations, we obtain for the 
expectation value of the trace (due to character orthogonality) 

/ i \ N—l °° 

RM^{Tr——)=-^^(-irv^e-^ c ^d(p, q ), (3.3) 

\ V W I p=Qq=0 

where C(p,q) and d(p,q) denote the value of the quadratic Casimir and the dimension of the 
irreducible representation identified by the Young diagram 





1 


2 






q 


1 

p 





C{p,q) = {p + q+\)(N- P + q N +l +q-p 
(N-p)(N + q) 1 



d(p,q)=d A (p)d\q) 



d\p) 



N p + q+l 1 



d S {q) 



We can exactly calculate sums of the form (with |z| < 1) 

N-\ 



£ u p d A {p){N-p) =N(l + i 

p=Q 



\N-\ 



^Z q d S {q){N + q) 

q=0 



N 



N+l " 



(3.4) 

(3.5) 
(3.6) 
(3.7) 

(3.8) 
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To factorize the sums over p and q in (3.3), we first write 

1 C l 



p + q + 1 Jo 



dpp p+q . (3.9) 



The ?-dependent weight factor is the exponent of a bilinear form in p and q (given by (3.5)). By a 
Hubbard-Stratonovich transformation the dependence on p and q can be made linear. Performing 
the (independent) sums over p, q then leads to [5] 

— , > N 2 , f r dxdy f l , w/JiJuif.,^ i (r fv Jl - vpe~ x ~'/ 2 ] N 1 

R(v) = e~2 / / / dp e -B(^+r)+s(*+w ^ L_ (3.10) 

f j i-» 271 7o [l-vpe'W 2 ] + 

(valid for |v| < 1). Now the entire dependence on N is explicit. The infmite-A^ limit of p^ ue can be 
obtained from this integral representation by using a saddle-point approximation for the integrals 
over x and y (cf. Sec. 5). 



4. Saddle-point analysis for p sym 

An integral representation for \jf(z) = (det(z — (which determines p sym ) is obtained in 

a similar manner [4]. In this case, only a single integral is needed (valid for \z\ > 1), 



NT N N 1 ( ;w ;»\- N 

V A(z) = e TJ__ J _Jwe-7r w (ze- l 2-e l 2) . (4.1) 
We set z = e' e+e and take the limit e — > + at the end. The integrand is exp(— Nf(w)) with 

2 

f{w) = ^+\og(ze-^-e^) . (4.2) 

Singularities of the integrand are located on the line Imw = — e < 0, which means that the integra- 
tion path (along Imw = 0) can be shifted upwards in the complex plane. 
The saddle point equation, f'(w s ) = 0, can be written as 

TU(e,T) U{e,T) + l/2 +i6 

l/(0,r)- 1/2 K ' 

with w s = iTU(6,T) = iT(U r (6,T) + iUi(6,T)). Taking the absolute value leads to 

Uf = U r coth{TU r + e) - U? - - . (4.4) 

This equation describes one or more curves in the complex-?/ plane (examples are shown in Fig. 2) 
on which the saddle points have to lie (it has been previously investigated in [6] for e = 0). For a 
given value of 6, the saddles are isolated points on these curves. 

For each d there is always one (and only one) saddle point on the closed curve encircling U = \ 
(cf. Fig. 3). The integration contour (along the imaginary-?/ axis) can be smoothly deformed to go 
through this saddle along a path of steepest descent (no singularities are crossed). 
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Figure 2: Examples of the contours in the complex-?/ plane described by equation (4.4) for T = 3 (left), 
T = 4 (middle), and T = 5 (right). The red dashed curves are for small e > 0, while the black solid curves 
are for e = 0. For our saddle-point analysis we keep £ > 0. 



Figure 3: Example for the location of the saddle 
points and the deformation of the integration contour 
in the complex-w plane for T = 5 and = 3 (with 
small e > 0). The thin solid lines are lines of con- 
stant Ref(w), the arrows point in the direction of in- 
creasing Re/(w). On each of the closed orange curves 
there is one saddle point (red dot and blue dot), and on 
the open orange curve there are infinitely many sad- 
dle points, but only one of them in the region shown 
in the plot (green dot). The thick blue curve is the in- 
tegration path along the direction of steepest descent 
through the relevant saddle point. 




The Gaussian integral which results from expanding the exponent f(w) to quadratic order 
around vt\ leads to 



U[ 1 + 



1 T(\/4-U 2 ) 
N[l-T{\/4-U 2 )] 2 



(4.5) 



where U = U(d, T) is the saddle point located on the closed curve around U = 1/2. Higher-order 
terms in the 1 /N expansion can be obtained by expanding the exponent to higher powers in w — w s . 
The leading-order result is 



limp^ ym (0,r) = 2Re[/(0,r), (4.6) 

which equals p oo (0, T) of Durhuus and Olesen [1] since U(d,T) satisfies (4.3). 

For T > 4, we find that Re U(d, T) is always non-zero (cf. right plot in Fig. 2). ReU(d,T) 
is also non-zero for T < 4 and \ 6\ < C {T). For T < 4 and \ 6\ = d c (T) we have f"(w s ) = 0, i.e., 
the first-order correction diverges and the Gaussian approximation does not work. For T < 4 and 
|0| > 6 C (T) the saddle point U (6,T) is purely imaginary (cf. left plot in Fig. 2) so that both the 
leading order and the 1 /N term are zero. In the region 1 6\ > 6 C (T) p^ ym is exponentially suppressed 
in ,/V and the study of the large-N asymptotic behavior requires more work. 

In Fig. 4 we show examples for the 1 /N corrections to Poo(0, T) for N =\0 and T = 2 and 5. 



5 



Infinite-N limit of the eigenvalue density of Wilson loops in 2D SU(N) YM 



Robert Lohmayer 





0.5 1.0 1.5 2.0 2.5 



Figure 4: Examples for the 1 /N corrections to p»o(0, T) for N = 10, T = 2 (left), and T = 5 (right). Shown 
are the exact result for p^ m (0,T) (blue dashed curve), the infinite-A^ result p 00 (9,T) (red dotted curve), 
and the asymptotic expansion of p^ m (9,T) up to order &{\JN) (black solid curve). We observe that the 
asymptotic expansion converges rapidly for small 1 9 1 and more slowly for larger 1 9 \ . 



true 



5. Saddle-point analysis for p 

The infinite-A^ limit of the true eigenvalue density p true can be obtained from a saddle-point 
approximation of the integral (3.10), which can be written in the form 



*(v) 



N- 



-e 2 



dxdy 



d p e -%{£+?)+h{x+iy) 2 -\(*-iy) 



t J J-oo 2% JQ 

x e {N-l) log( l-vpe-*->/ 2 ) - (N+ 1 ) log( 1 -vpe-y-'l 2 ) 
i0-e 



(5.1) 



(valid for |v| < 1; we set v = e and take the limit e — > + after deforming the integration 
contour). We first approximate the integrals over x and y and integrate over p at the end. At 
leading order, the integrals (over x and y) decouple and can be approximated independently. The 
coefficients of — N in the exponent are 



f(y) 



1 

2r 
1 

2r 



y + log 

log 



1 - vpe'- 
1 — vpe 



-x-t/Z 



-f{ix) 



(5.2) 
(5.3) 



Let us first consider the integral over y. Substituting y = w — it /2 leads to the integrand for y 
in (4.1) (with z — ► 1/vp). The integration over w is now along the line from —°° + it/2 to +oo + #/2 
(since there are no singularities between this line and the real-w axis, the integration contour can 
be shifted to the real-w axis). The saddle-point equation reads (with y s = it(U — 1/2)) 

v U+l/2 = 1 
U-l/2 vp ' 



(5.4) 



The only difference to Eq. (4.3) is that the absolute value of the RHS can take any value in the 
interval < \vp\ < 1 (for p = 1, we are in the situation of Sec. 4). Numerical analysis shows that 
the integral can again be approximated by one single saddle point, yo(6,t,p). The relevant saddle 
point is again located on a closed curve around [7 = 1/2 (corresponding to y = 0). This curve 
shrinks with decreasing p (cf. Fig. 5). 

Due to / (x) = — / (be) , the saddle points of the x- and j-integrals are related by a rotation of 
7i/2 in the complex plane, x s = —iy s . The relation to U is x s = t(U — 1 /2), the integration path for 
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Figure 5: Contours of solutions of equation (5.4) in the complex-C/ plane at t = 3 (left), t = 4 (middle), and 
t = 5 (right) for p = 1 (black, solid), p = 0.9 (red, dashed), p = 0.6 (green, dot-dashed), and p = 0.3 (blue, 
dotted). 



the .x-integral is along the real-£/ axis. By analyzing the directions along which the phase of the 
integrand is constant, we find that the integration contour can always be deformed to go through 
the (single) saddle point in the right half-plane (on the curve around U = 1/2) in the direction of 
steepest descent. Depending on p, v, and t, there is either one or no additional saddle point on the 
contour(s) in the left half-plane through which we can also go in the direction of steepest descent, 
but the contribution of this additional saddle point (if there is one) is exponentially suppressed. 
Therefore the relevant saddle point is xq = —iyo, i.e., both integrals can be approximated by a 
single saddle point. An example for the location of the saddle points and the deformation of the 
integration path is shown in Fig. 6. 



Figure 6: Example for the location of the saddle 
points and the deformation of the integration path in 




the complex-t/ plane for t = 5 and p = 0.95. The 
dashed black curves (two closed, one open) are the 
curves on which all saddle points have to lie. In this 
example 6 = 3.0. On each of the closed curves there 
is one saddle point (red dot and blue dot), and on 
the open curve there are infinitely many saddle points, 
but only one of them in the region shown in the plot 
(green dot). The thin solid lines are lines of constant 
Re/(x) and Re/(y). The arrows point in the direc- 
tion of increasing Re/(x) or decreasing Re/(y). The 
dashed blue curve is the integration path for the y- 



integral along the direction of steepest descent. The 
solid red-blue curve is the integration path for the x- 
integral along the direction of steepest descent. 



Combining the saddle-point approximations for the x- and y-integrals leads to (xq =xo(6,t,p)) 

t Jo t+x (t+x ) 

By differentiating the saddle-point equation with respect to p we find that the integrand is just 
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v 1 -j^ and therefore 

N T y dxn N 
R(v) = -- dp^ = --[xo(e,t,p = l)-x (e,t,p = 0)}. (5.6) 
tv Jo dp tv 

The eigenvalue density at infinite N is determined by the locations of the saddle points at p = 1 
and p = only. Since xo(d,t,p = 0) = 0, we obtain 

limp7 e (M =2ReU(6,t,p = 1) = lim p s N ym (0J), (5-7) 

which confirms that the infinite-Af limit of the true eigenvalue density is given by p cx ,(d,t) of 
Durhuus and Olesen [1]. 

6. Conclusions 

The probability distribution of Wilson loops in SU(A^) YM in two Euclidean dimensions can 
be written as a sum over irreducible representations (where only dimension, second-order Casimir, 
and character of W enter). This allows for the derivation of integral representations for different 
density functions (including the true eigenvalue density), which have the same infinite-A^ limit. 
These integral representations, where N enters only as a parameter, are exact for any finite N. 
Results at infinite N can be obtained by saddle-point approximations. Next-order terms in an 
expansion in 1 /N give reasonable results in the interval where p^Q) > 0. More work is needed to 
get finite-A^ effects where p«,(0) = and corrections are exponentially suppressed in N. 
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